Theory of giant and nil proximity effects in cuprate semiconductors 
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A number of observations point to the possibility that high-T c cuprate superconductors may not 
be conventional Bardeen-Cooper-SchriefTer (BCS) superconductors, but rather derive from the Bose- 
Einstein condensation (BEC) of real-space pairs, which are mobile small bipolarons. A solution of 
the Gross-Pitaevskii equation describing bose-condensate tunnelling into a cuprate semiconductor 
is analytically found. It accounts qualitatively and quantitatively for nil and giant proximity effects 
discovered experimentally in cuprates. 
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Perhaps the most striking property of cuprate super- 
conductors is the Giant Proximity Effect (GPE). Sev- 
eral groups Q reported that in the Josephson cuprate 
SNS junctions super current can run through normal TV- 
barriers with the thickness 2L greatly exceeding the co- 
herence length, when the barrier is made from a slightly 
doped non-superconducting cuprate (the so-called N' 
barrier). Using an advanced molecular beam epitaxy, 
Bozovic et al 1] proved that GPE is intrinsic, rather 
than caused by any extrinsic inhomogeneity of the bar- 
rier. Resonant scattering of soft-x-ray radiation did not 
find any signs of intrinsic inhomogeneity (such as charge 
stripes, charge-density waves, etc.) either Q. Hence 
GPE defies the conventional explanation, which predicts 
that the critical current should exponentially decay with 
the characteristic length of about the coherence length, 
£ < 1 nm in cuprates. Annealing the junctions at low 
temperatures in vacuum rendered the barrier insulating. 
Remarkably when the SN'S junction was converted into 
a superconductor-insulator-superconductor (SIS) device 
no supercurrent was observed, even in devices with the 
thinnest (one unit cell thick) barriers 3( (nil proximity 
effect, NPE). 

Cuprate superconductors with low density of free car- 
riers and poor mobility (at least in one direction) are 
characterized by poor screening of high-frequency c-axis 
polarised optical phonons. The unscreened Frohlich in- 
teraction between oxygen holes and these phonons com- 
bined with on-site repulsive correlations (Hubbard U) 
binds holes into superlight intersite bipolarons 4], which 
are real-space pairs dressed by phonons. Experimental 
evidence for exceptionallystrong e-ph interactions is now 
so overwhelming [Mjp, 0, 0, LJ] that a bipolaronic charged 
bose gas (CBG) [l(| could be a feasible alternative to 
the BCS-like scenarios of cuprates. The bipolaron the- 
ory predicted such key features of cuprate superconduc- 
tors as anomalous upper critical fields, spin and charge 
pseudogaps, and unusual isotope effects later discovered 
experimentally. The theory explained normal state kinet- 
ics, including the anomalous Hall-Lorenz number, high T c 
values, specific heat anomalies of cuprates (for a review 
see Ul|), and more recently the d-wave symmetry of the 



order parameter in underdoped and overdoped [T^ 
samples, the normal state Nernst effect 13] and diamag- 
netism |l4j . 

Here I show that both GPE and NPE can be broadly 
understood as the BEC tunnelling into a cuprate semi- 
conductor. 

A stationary condensate wave function of CBG, ^(r), 
obeys the following equation, 
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^(r) = 0, 



(1) 



where m is the boson mass, \i is the chemical potential 
and V(r) describes the Coulomb and hard-core composed 
boson-boson repulsions, and their attraction to a neutral- 
ising charge background (here and further H = c = ks = 
1). While the electric field potential can be found from 
the corresponding Poisson-like equation ^| , a solution of 
two coupled nonlinear differential equations for the order 
parameter ^(r) and for the potential V(r) is a nontrivial 
mathematical problem, which requires a " multishooting" 
numerical approach. Here we restrict our analysis by 
a short-range potential, V(r) = V^r)! 2 , where a con- 
stant V accounts for the hard-core and screened Coulomb 
repulsions. Then Eq.(l) is the familiar Gross-Pitevskii 
(GP) equation [l6|. In the tunnelling geometry of SN'S 
junctions, Figs. 1,2, it takes the form, 



1 d 2 ^(Z) 
2m~ c dZ 2 



[V\il>(Z)\ 2 - M(Z), 



(2) 



in the superconducting region, Z < 0, near the SN 
boundary, Fig.l. Here m c is the boson mass in the direc- 
tion of tunnelling along Z. Deep inside the superconduc- 
tor \ip(Z)\ 2 = n s and \i = Vn s , where the condensate 
density n s is about x/2, if the temperature is well below 
T c of the superconducting electrode. Here the in-plane 
lattice constant a and the unit cell volume are taken as 
unity, and x is the doping level as in I^-xSr^CuO^ 

The normal barrier at Z > is an underdoped cuprate 
above its transition temperature, T' c < T where the 
chemical potential \i lies below the bosonic band by some 
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FIG. 1: BEC order parameter at the SN boundary for 
fl = 1.0,0.1,0.01 and ^ 0.001 (upper curve). The chemi- 
cal potential is found above the boson band-edge due to the 
boson-boson repulsion in cuprate superconductors and below 
the edge in cuprate semiconductors with low doping. 



energy e, Fig.l, found from 

J dEN(E)[exp((E + e)/T) - l]" 1 = x f /2. (3) 

Here N(E) is the bipolaron density of states (DOS), and 
x' < x is the doping level of the barrier. In-plane bipo- 
larons are quasi- two dimensional bosons, propagating 
along the CuO^ planes with the effective mass m several 
orders of magnitude smaller than their out-of-plane mass, 
m c ^> m, which accounts for the remarkable anysotropy 
of cuprates at low temperatures H?|- Using bipolaron 
band dispersion, E(K) = K 2 /2m + 2t c [l — cos(K±d)], the 
density of states is found as N(E) = (ra/27r 2 ) arccos(l — 
E/2t c ) for < E < 4t c , and N(E) = m/2n for 4t c < E. 
Here K and K± are the in-plane and out-of-plane center- 
of-mass momenta, respectively, t c = l/2m c d 2 , and d 
is the inter-plane distance. Assuming that the critical 
temperature of the barrier T' c « To/ ln(To/2£ c ) (here 
To = irx f /m ^> T' c ) is much larger than t c one obtains 
from Eq.(3) 

e(T)^-Tln(l-e- T °/ T ), (4) 

which turns into zero at T = T' c . 

Then the GP equation in the barrier is written as 

1 d 2 i){Z) 



[VmZ)\ 2 + eMZ). 



(5) 



2m c dZ 2 

Introducing the bulk coherence length, £ = 
l/(2m c n s V) 1 / 2 and dimensionless f(z) = ip(Z)/nl^ 2 , 
fl = e/n s V, and z = Z/£ one obtains for a real /(z) 

fl 

dz 2 



f ~ f, 



(6) 



(7) 



if z > 0. These equations can be readily solved using first 
integrals of motion respecting the boundary conditions, 
/(-oo) = 1, and /(oo) = 0, 



dz 



= -{i/2+t/2-n 



2U/2 



and 



dz 



\l/2 



(8) 



(9) 



for z < and z > 0, respectively. The solution in the 
superconducting electrode is given by 



f(z) = tanh 



-2 1/2 z + 0.51n 



2 1 ' 2 (l + p,) 



1 



2 1 /2(l + /i )l/2_ 1 



(10) 



It decays in the close vicinity of the barrier from 1 to 
/(0) = [2(1 + /i)] -1 / 2 in the interval about the coherence 
length £. On the other side of the boundary, z > 0, it is 
given by 



/(*) 



(2£) 



sinh^/a + ln[2(A(l + /2)) 1 / 2 + (1 + 4/i(l + /i)) 1 ^]} ' 

(11) 

Its profile is shown in Fig.l. Remarkably, the order pa- 
rameter penetrates the normal layer up to the length 
Z* « (/i) _1 / 2 £, which could be larger than £ by many 
orders of magnitude, if fl is small. It is indeed the case, if 
the barrier layer is sufficiently doped. For example, tak- 
ing x' = 0.1, c-axis m c = 2000m e , in-plane m = 10m e 
[III, a = 0.4 nm, and f = 0.6 nm, yields T « 140 K and 
(/i) -1 / 2 > 50 at T = 25K. Hence the order parameter 
could penetrate the normal cuprate semiconductor up to 
a hundred coherence lengths or even more as observed 
Q. If the thickness of the barrier L is small compared 
with Z*, and (fl) 1 / 2 <C 1, the order parameter decays 
following the power law, rather than exponentially, 



/(*) 



V2 
z + 2' 



(12) 



Hence, for L < Z*, the critical current should also de- 
cay following the power law as discussed below. On the 
other hand, for an undoped barrier fl becomes larger than 
unity, fl oc ln(raT '/ttx') — > oo for any finite temperature 
T when x' — > 0, and the current should exponentially 
decay with the characteristic length smaller that £, as is 
experimentally observed as well 0. 

To get more insight in the temperature and barrier- 
thickness dependence of the critical current one has to 
solve the GP equation in the SN' 'S junction geometry, 
Fig. 2, with the current, J = (2en s /m c ^)R 2 dQ/dz. Here 
R(z) and O(z) are the amplitude and the phase of the 
order parameter, respectively, /(z) = R(z) exp(i6(z)). 
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FIG. 2: Numerical values of the intgral I(p) in Eq.(19) (sym- 
bols) compared with the analytical approximation (solid line) . 
Inset shows schematically the trilayer SN f S device of Ref.[l|. 



The current does not depend on z, so that the equation 
can be written as 



d 2 R 
~dz^ 



if \z\ > Z, and 



d 2 R 
~dz^ 



R 3 -(l-f/R*)R, 



R 3 + (p + j 2 /R 4 )R, 



(13) 



(14) 



if \z\ < Z, where Z = L/£ and j — Jm c £ ) /2en s . Using 
integrals of motion respecting the boundary conditions, 
dR/dz = at z = ( R(z) = R{—z) since the junction 
is symmetric), R 3 ^ = R^ — j 2 / R^, and dR^jdz = 0, 
these equations are reduced to 



dR 

dz 

for < z < Z, and 
dR 



[(R* - Rt)/2 + f,(R 2 - R 2 ) + j 2 {Rf - ir 2 )] 1 /*, 

(15) 



dz 



l(R 4 - RD/2 + (Rl - i? 2 ) + j\R~J - i?" 2 )] 1 / 2 , 

(16) 

for z > Z. Then integrating Eq.(14) from z = up to 
z = Z allows us to connect the current and the order 
parameter, Rq in the center of the barrier at z = 0, 



dx 



y / x(x 2 - 1) + xjx - l)a +(x- 1)13 



V2R L 
(17) 

Since Ri ^> Ro for sufficiently thick barrier, Z ^> 1, the 
upper limit in this integral can be extended to oo yielding 
a transcendental equation for Rq as a function of j 



2^1 



1/2, 1/2; 1; 1(1 - ^-J=^=^)' 



V / 2(a + / 9 + 2) 1 / 4 i?o«A. 



(18) 



where 2Fi(a,b;c; z) is the Gauss hypergeometric func- 
tion, a = and /? = 2j 2 /R$. When a = /? = 0, one 
obtains Ro oc 1/Z, as in Eq.(12), and when a — > oo, the 
order parameter decays exponentially as in Eq.(ll). 

To evaluate the critical current one can reduce the in- 
tegral in Eq.(17) to a single parameter integral, I(p) = 
J °° dx[x(x^l) 2 ^px}-^ 2 withp= -l+4(a+/3+2)/(3 + 
a) 2 . This integral can be analytically approximated (see 
Fig.2) as I(p) « 3.3(1.1 +p)" 1/4 , if p is not too close to 
— 1. Then solving 



I(p) =R l(3 + a) 



1/2 



(19) 



with respect to the current yields j ~ Ro(25/2l 4 — Rq — 
pRl) 1 / 2 , and the critical current j c « 4(1 - 2pl 2 )/l 3 . 
These expressions are applied when I < (/i) -1 / 2 (i.e. L < 
Z*). 

To get the temperature dependence of the critical cur- 
rent J C (T) in the whole range of parameters one can ap- 
ply the scaling, j c oc Rq/1, where Rq is defined as the 
exact solution of the SN' boundary problem, Eq.(10) at 
z = Z, Rq = f(l). In such a way one finds (in ordinary 
units) 



Jc{T) 



Ae(T) 



k B smYi 2 [L^2m c e{T)/h] ' 



(20) 



where A w AekB^QTiso/hL is temperature independent. 
With the zero-temperature coherence length £o = 0.5 
nm and condensate density n s o = x/2Q one estimates 
A ^ 150 kA/cm 2 K for optimum doping x = 0.15 and 
the barrier thicness 2L = 10 nm (here Q is the unit 
cell volume). Eq.(20) fits well the temperature depen- 
dence and absolute values of the critical current, Fig. 3, 
measured by Bozovic et al. [l| with To = 160 K, and 
m c = 1000m e , Fig. 3. The typical device of Ref.0 used 
Lai.85Sro.i5Cu04 with T c w 45K as the superconducting 
electrodes while the normal N f barrier was made of un- 
derdoped La2Cu04+5 with typical T' c w 25K. One can see 
from Fig. 3 that the theory describes the critical current 
in the normal region of the barrier, where jl is positive. 
When the barrier becomes superconducting (c.a. below 
20 K) the experimental critical current naturally exceeds 
the theoretical "semiconducting" J C (T). 

A possibility of real-space pairing, proposed originally 
by Qgg 1 181 and later on by Schafroth and Blatt and 
Butler |l3 has been the subject of many discussions 
as opposed to the Cooper pairing, particularly heated 
over the last 20 years after the discovery of high tem- 
perature superconductivity in cuprates. Our extension 
of the BCS theory towards the strong interaction be- 
tween electrons and ion vibrations proved that BCS and 
Ogg-Schafroth pictures are two extreme limits of the 
same problem. For a very strong electron-phonon cou- 
pling, polarons become self-trapped on a single lattice 
site and bipolarons are on-site singlets. In the Holstein 
model of the electron-phonon interaction their mass ap- 
pears only in the second order of polaron hopping, so 
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FIG. 3: Temperature dependence of the critical current of the 
trilayer SN'S device 1] with L = 5 nm (symbols) described 
by Eq.(20) at T > T' c . 



that on-site bipolarons are very heavy. This estimate 
led some authors to the conclusion that the formation 
of itinerant small polarons and bipolarons in real mate- 
rials is unlikely |20|, and high-temperature bipolaronic 
superconductivity is impossible [2l[. However, we have 
noticed that the Holstein model is an extreme polaron 
model, and typically yields the highest possible value 
of the (bi)polaron mass in the strong coupling regime. 
Cuprates are characterized by poor screening of high- 
frequency optical phonons and are more appropriately 
described by the long-range Frohlich electron-phonon in- 
teraction |4j. The unscreened Frohlich electron-phonon 



interaction provides relatively light small polarons and 
bipolarons, which are several orders of magnitude lighter 
than small Holstein (bi)polarons. 

I conclude that the bipolaron theory accounts for GPE 
and NPE in slightly doped semiconducting and undoped 
insulating cuprates, respectively. It predict s the occur- 
rence of a new length scale, h/ ^/2m c e(T). In a wide 
temperature range far from the transition point, T' c < 
T < Tq, this length turns out much larger than the zero- 
temperature coherence length, if bosons are almost two- 
dimensional (2D). The physical reason, why the quasi- 
2D bosons display a large normal-state coherence length, 
whereas 3D Bose-systems (or any-D Fermi-systems) at 
the same values of parameters do not, originates in the 
large DOS near the band edge of two-dimensional bosons. 
Since DOS is large, the chemical potential is pinned near 
the edge with the magnitude, e(T), which is exponentially 
small when T < Tq. Importantly the theory predicts an 
unusual dependence of J c on the barrier thickness 2L, in 
particular at low temperatures, T < T' c1 when the bar- 
rier is superconducting, there is almost no L dependence, 
above the transition of the barrier, J c oc 1/L 3 , and well 
above the transition (T > To >> T' c ) J c decays exponen- 
tially with L. Based on a great number of experimental 
observations 10] including GPE and NPE I argue that 
the most likely scenario for superconducting cuprates is 
the genuin Bose-Einstein condensation of real-space mo- 
bile lattice bipolarons. 
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